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Abstract 

We present a model of baryonic matter which contains free con- 
stituent quarks in addition to bound constituent quarks in nucleons. 
In addition to the common linear cj-model we include the exchange of 
vector-mesons. The percentage of free quarks increases with baryon 
density but the nucleons resist a restoration of chiral symmetry. 
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1 Introduction 



The possibility of creating high density baryonic matter has triggered several 
speculations about the appearance of a quark phase and/or diquark phase 
with or without chiral symmetry restoration. Most naive quark models with 
effective interactions have an inherent very low scale at which chiral sym- 
metry is restored. This scale can be obtained from an effective potential 
calculation within a linear sigma model for example, and varies between 0.5 
and 1.5 times normal nuclear density po- This density appears to be unre- 
alistically low, because higher mass mesons (e.g. a;-mesons) are neglected in 
these models. The inclusion of higher mass mesons yields repulsive forces 
that drive the system away from the chiral symmetry restoration phase tran- 
sition. 

Relativistic Quantum Hadron Dynamics (QHD) has been a very successful 
mean field theory of nuclear matter for certain features of nuclei. It gives 
a description of nuclear saturation and of the spin-orbit potential in nuclei, 
which is close to that obtained phenomenologically. On the other hand, its 
compressibility K^o > 500 MeV is too high, compared to the empirical value 
= 210 MeV [|l|. Although due to confinement nucleon degrees of freedom 
replace quark degrees of freedom in vacuum, a growing part of the nucleon 
wave functions will overlap with increasing nuclear density, because of the 
finite size of the nucleons. By this, the substructure of the nucleon matters 
already before any phase transition, and so it is interesting to ask the ques- 
tion what chemistry of nuclear matter results if one combines both nucleons 
and quarks. How is baryonic matter built up when a fraction of baryonic 
charge can be delocalized in quarks, which have their own dynamics? In the 
seventies various schemes of covalent bonding were invented. The above idea 
goes far beyond these. It assumes the existence of an unconstrained quark 
substrate. 

2 Quark nucleon model 

We unify both quark and nucleon dynamics by a common linear cr- mo del 
Lagrangian via which both sectors communicate. In equilibrium we do not 
need complicated transition terms n ^ 3q. 
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The Lagrangian of our model reads 

+ \ {{D'^r^^r + (Dl^^'^cTr) - ^n.^'^-^ - via' + n^) , 

with the covariant derivatives D^i'"'^'^'"'^ = d^ + igyn/sv^ii, the field tensor = 
d^ujy — dyUJ^ and the spontaneous symmetry breaking potential (m^ < 0) 

V(cT +7r) = — (a +7r) + -(a + tt ) . 

The free constituent quarks with the same mass as quarks bound in nucleons 
represent larger subclusters of quarks, for which a definite baryonic center 
point can no longer be localized. These "free" quarks mimic the frequent 
switching of string configurations. In this picture, the a;-repulsion is identified 
with the repulsion of the string junctions. 

The hadronic part is a modification of the model presented by Walecka [0], 
with the explicit mass term of the vector meson replaced by a coupling to 
the scalar mesons. Thus in mean field theory, the uj also gains its mass from 
the (T-field. This leads to a scaling of all hadron masses as suggested by 
Brown and Rho [Q], who have also investigated the compatibility of nucleon 
and quark degrees of freedom at finite temperature and density [Q. As a 
consequence, the chiral symmetric phase with (a) ~ 0, i.e. very small u- 
mass, is suppressed, because in the mean field approximation the a;-repulsion 
becomes very large. 

We choose the following couplings: ggq = 3.23, gsn = ^ ■ Qsq = 9.69, ml = 
—0.26, A = 30.0, gsy = 8.41, g^n = 9.5. The couplings of the mesonic 
potential and A are calculated from a renormalization group (RG) fiow- 
equation approach for the quark part of the model 0. The couplings yield 
a cr- mass of nicr = 720 MeV and a nucleon mass of m„ = 901 MeV. The 
nucleon mass is less than 938 MeV because an explicit symmetry breaking 
term is not included in our Lagrangian. In the mean field approximation 
with constant expectation values a and Co, the Lagrangian takes the form 

^MF = i^ni^l^d^ - g^nl^U) - gsn^)4'n + IpqilY - gsq^)i^q 
I 1 2-2-2 '^s -2 ^-4 

+ ^9s.cru -—a --a . 
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The number of internal degrees of freedom is 7„ = 4 for the nucleons and 7^ 
= 12 for the quarks, because of 2 spin and 2 isospin states for both and in 
addition 3 different colors for the quarks. In the ground state, the nucleons 
and quarks build two Fermi gases with Fermi momenta kpn and kpg. They 
are connected to the conserved overall baryon density Pb = 'ipni'n + ^'4'l''Pg 
through 



x)pB = r^^Fn and 3xpB = , 



with X denoting the concentration of quarks in the system. 
Now one can calculate the energy density, which isQ 




kin + 9L^' 4n + ^^^^^V^ 1 - ^ log 



9sn'^^kFn\ gtn(^^ , „ ( kpn + V kpn + dsn^^^ 



+ ^ {kFn - k,„ ^ + T^f-^^n + + , 

where the mean field value a in the ground state is self-consistently given by 
minimizing e with respect to a. 



3 Results 

The result of our calculation is plotted in Fig.|I|. To the left the x = 
boundary curve shows the pure hadronic phase. The right boundary at a; = 1 
shows the same for pure quark matter with a kink due to the first order chiral 
phase transition. This transition persists only for a quark concentration x > 
0.99, where at the endpoint the phase transition is of second order. The white 
line denotes the stable configuration with minimal energy at each density 
and shows a smooth transition between a vanishing quark concentration in 
vacuum and an asymptotic quark concentration x^o ~ 0.9 as pb oc. 
This curve of the energy minimum, which represents the equation of state 
(EOS) of nuclear matter, is also plotted as a 2D projection in Fig.|| together 
with the advanced nuclear physics calculations by Pandharipande et. al. [^]. 
Whereas quark matter and nuclear matter alone do not bind, the mixture of 
both binds at a saturation density ps = 1.14po with Eb = 16 MeV. 

^The first term in the second line is the same as the nucleon term in the first line with 
interchanged momenta and couplings 
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Figure 1: Binding energy per baryon ^ = "^^^"^j"^ — rriN as a function of 
the independent variables x and ps given in units of normal nuclear matter 
density po = 0.17 fm~^ 

The Dirac mass of the nucleon at saturation is m*^ = Qsn^^s = 581 MeV. 
As pointed out in this quantity is not the relativistic analogon of the 
nonrelativistic effective mass measured in nuclear experiments. The corre- 
sponding quantity is the relativistically equivalent effective mass m*. We get 
^ = 1 - ^ = 0.71 compared to an experimental value oi — ^ 0.74 fZH. 
As in the standard Walecka model the compression modulus K^o = 823 MeV 
is much too large. This is a general problem of these simple mean field mod- 
els and could be corrected by adding further meson interactions or potential 
terms to the Lagrangian. 

The saturation value for the quark concentration is Xg = 50%. This rather 
high value is reduced in finite nuclei due to the large surface of low nuclear 
matter density. A simple estimation with a Fermi distribution for the density 
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Figure 2: Nuclear matter equation of state 



gives (x) ~ 40% for ^osp]^ g^j-^^ ^2;) < 30% for ^^O. There are indeed exper- 
iments on the occupancies of lower shell model orbitals, which suggest that 
these orbitals are not entirely filled. Instead the Fermi surface is smeared 
out. This effect has been explained by a mixing of shell model orbitals, but 
in our model it is simply due to the coupling of different quark momenta to 
a baryonic three quark system. Such a composite baryon is not necessarily 
a nucleon but may be another excited baryonic state. Fig.| shows the occu- 
pancies in lead computed with this model and looks almost the same as the 
result obtained under the inclusion of ground state correlations [Q. In the 
model described here, the "free" constituent quarks reduce the amount of 
quasi baryons below the Fermi momentum and enhance their number with 
momenta higher than kpn- The unchanged part of the Fermi surface is con- 
nected to the residue factor z, which can be measured. Because the 3s-orbital 
in lead is nearly exclusively located in the center of the nucleus, one can take 
in very good approximation the nuclear matter value of the quark concen- 
tration to compute the residue factor for this orbital. The result of 2; = 0.67 
compares favorably to the value z = 0.64 ± 0.06 from (e, e'p) experiments 
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Figure 3: Occupation numbers n(k) of quasi baryons in lead 

4 Conclusion 



Because of the repulsive cj-mean field at saturation gynCOg = 263 MeV, the nu- 
cleons occupy the higher energy levels, whereas the constituent quarks have 
a strong occupancy in the deep lying levels. In such a way the overlap of con- 
stituent quarks and nucleons is restricted to the upper levels and would lead 
to a small correction of the energy density. Note that this overlap remains 
small at densities regarded here. The validity of our calculation is also lim- 
ited by the credibility of the mean field approximation and the appearance 
of gluonic degrees of freedom at high density. We estimate the momentum 
scale for this region to be fc[/y > 1 GeV cf. with our RG calculation [|1^ 



Due to the nonabelian nature of QCD the three body gluon exchange forces 
are more important than the two body gluon exchanges. The three body 
forces lead to baryonic tripling whereas two body forces induce diquark for- 
mation |Tl]]. The latter may be an intermediate stage on the way from high 
density to low density, but as shown in nature color neutralization will win 
ultimately at lower densities. 
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To summarize, we have constructed a very simplified model with constituent 
quarks and nucleons, which interpolates between low density nucleon matter 
and high density quark matter in a continuous way. As a result we have 
shown that our model qualitatively reproduces the nuclear matter EOS and 
provides a simple explanation for experimental measured residue factors. 
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